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基于 小 波 变换 模 的 FAB 扩散 模型 
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(西安 交通 大 学 理学 院 信 息 与 计算 科学 系 ， 西 安 710049) 

摘 E 将 小 波多 尺度 边缘 检测 的 思想 引入 到 前 向 一 后 向 (FAB) 扩散 模型 中 ， 根 据 小 波 变换 能 够 准确 提 
取 图 像 局 部 结构 信息 的 特性 ， 给 出 了 利用 多 个 尺度 上 小 波 变 换 模 构造 图 像 扩 散 系 数 的 方法 ， 进 而 
建立 了 基于 小 波 变换 模 的 前 向 一 后 向 扩散 模型 . 最 后 给 出 了 一 系列 数值 实验 结果 ， 以 比较 本 文 方 
法 与 传统 FAB 方法 的 图 像 平滑 增强 效果 . 
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近年 来 ， 基 于 偏 微分 方程 (PDE) KRAAL] 受到 了 广泛 的 关注 . 早期 的 工作 主要 包括 
由 Witkin 引入 的 线性 尺度 空间 理论 叫 ， 该 理论 将 扩散 过 程 的 扩散 时 间 与 高 斯 平滑 滤波 器 的 方 
差 联 系 起 来 ， 建 立 了 线性 扩散 方程 与 线性 滤波 系统 的 等 价 性 . 1990 年 ，Perona 和 Malik 引入 了 
非 线 性 扩散 方法 ， 常 称 为 P-M 扩散 模型 中 ,该 模型 的 扩散 系数 是 关于 梯度 模 的 非 负 下 降 函 数 ， 
它 使 扩散 沿 着 梯度 的 正方 向 进行 ， 并 且 使 扩散 速度 随 着 梯度 模 的 增 大 而 减 小 ， 由 于 平缓 区 域 的 
梯度 模 小 ， 扩 散 速度 大 ， 而 边缘 的 梯度 模 大 ， 扩 散 速 度 小 ， 因 而 P-M 模型 可 在 去 除 噪 声 的 同 
时 能 很 好 地 保持 边缘 . 但 P-M 模型 易 受 噪声 干扰 ， 而 且 是 不 适 定 的 . 1992 年 ，Catt6 基于 高 斯 
预 滤波 方法 改进 了 P-M 模型 由，Catt6 模 型 是 适 定 的 ， 而 且 由 于 采用 了 高 斯 梯度 ， 扩 散 系数 的 
估计 较为 准确 ， 从 而 提高 了 恢复 图 像 的 质量 . 然而 P-M 模型 和 CattE 模型 均 采 用 了 同样 性 质 的 
扩散 函数 ， 使 得 扩散 过 程 在 边缘 附近 会 减 小 或 停止 扩散 ， 导 致 边缘 附近 的 噪声 在 很 大 程度 上 被 
RA. 为 了 克服 这 一 缺陷 ，Weickert 引入 了 扩散 张 量 的 概念 ， 建 立 了 基于 扩散 张 量 的 平滑 增强 
AE], Weickert 模型 利用 了 梯度 方向 与 梯度 模 的 信息 ， 并 从 扩散 方向 与 扩散 量 两 方面 控制 扩 
散 过 程 ， 从 而 使 得 扩散 滤波 可 沿边 缘 进行 . 但 Weickert 方法 存在 过 度 扩 散 现 象 ， 常 使 恢复 图 像 
产生 奇异 波动 ， 故 常用 于 诸如 指纹 一 类 的 流 线 状 图 像 平滑 ， 而 不 用 于 自然 图 像 的 滤波 处 理 . 

上 述 扩散 模型 均 采 用 前 向 ( 正 向 ) 扩散 对 图 像 进行 处 理 ， 虽 然 对 边缘 有 保护 作用 ， 但 无 增强 
作用 . 为 了 在 去 噪 的 同时 增强 边缘 ，2002 年 ，Gilboa 提出 了 前 向 一 后 向 扩散 (FAB) RHE, 
该 模型 具有 平滑 增强 作用 ， 它 将 扩散 分 为 正 、 反 两 个 方向 的 扩散 : 正 向 扩散 起 平滑 作用 ， 反 向 
扩散 起 锐 化 边缘 的 作用 . 一 幅 图 像 经 前 向 -后 向 扩散 后 ， 灰 度 值 常 超出 初始 图 像 灰 度 值 范围 ， 
因而 不 满足 “ 极 值 原理 ”. 但 Gilboa 在 文献 [6] 中 指出 ， 在 局 部 范围 内 使 用 后 向 ( 闭 ) 扩散 不 会 
破坏 扩散 模型 整体 的 稳定 性 .我们 知道 ， 经 典 的 图 像 平滑 与 增强 是 独立 进行 的 两 种 处 理 过 程 ， 
而 FAB 扩 散 方法 则 将 二 者 结合 在 一 起 进行 ， 因 此 有 关 FAB 扩散 滤波 的 研究 在 近年 来 受到 了 人 
们 的 广泛 关注 [73]. 

本 文 的 研究 目的 是 基于 小 波 变 换 建 立 稳健 的 FAB 扩 散 模 型 . 一般 来 说 ， 扩 散 系 数 的 确定 
分 为 两 步 : 一 是 检测 边缘 信息 ， 二 是 确定 扩散 函数 . 其 中 边缘 信息 对 扩散 量 ( 即 扩散 系数 的 大 
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小 ) 的 确定 有 非常 重要 的 影响 . 在 Gilboa 所 给 的 FAB 扩散 模型 中 ， 边 缘 信 息 是 用 图 像 梯 度 检 测 
的 ， 这 使 得 扩散 过 程 对 噪声 比较 敏感 . 为 此 ， 本 文采 用 基于 小 波 变换 方法 个 提取 图 像 局 部 结构 
信息 ， 以 降低 噪声 对 扩散 系数 的 影响 ， 提 高 图 像 的 扩散 滤波 效果 ， 


2 预备 知识 
2.1 FAB 模型 
Gilboa 提出 的 FAB 扩散 模型 如 下 
| afeut) = div[C(\VF)VF], te (0,400), (z,y) €R, 
f(z,y,0) = folz, y), 
R (1) F fo(z,y) ASRA, Cls)\ 称 为 扩散 函数 ， 常 取 其 形式 如 下 


1 a 
TIFE TT ba 中 


其 中 ks, ko, w, a, n,m 为 参数. 参数 ky 的 作用 与 P-M 模 型 中 参数 作用 类 似 ， 用 来 定义 前 癌 
扩散 范围 ， 有 具体 见 文献 [3}， 参 数 如, w 用 来 定义 后 向 扩散 范围 ，ks 为 后 向 扩散 的 中 心 梯度 
E, wA k 为 中 心 的 后 向 扩散 宽度 ， 铝 的 取 值 取决 于 所 要 锐 化 增强 的 梯度 模 的 大 小 ， 并 且 通 
过 ww 来 控制 后 向 扩散 的 范围 ， 使 前 后 向 扩散 范围 不 会 发 生 重修， 参数 a 用 来 确定 前 后 向 扩散 比 
例 . 一 般 取 n = 4, m = 1， 而 其 它 参 数 应 满足 如 下 两 个 不 等 式 四 


C(s) 





ky 

< = 

wS a PaT 0<w< kp kf, (3) 
nee {s CH = Rieter tie is C(s)}; 4) 


AF s- C(s) 称 为 通 量 . 条 件 (4) 保证 在 后 向 扩散 过 程 中 不 会 产生 新 的 奇异 点 ， 从 而 不 破坏 扩散 
模型 的 整体 稳定 性 . 

FAB 扩 散 模 型 的 缺点 是 参数 很 难 确定 . Gilboa 通 过 计算 平均 梯度 绝对 值 (mean absolute 
gradient, MAG)， 即 |Vul 的 平均 值 来 确定 参数 ， 例 如 : [ky, kow] = [2,4,1] x MAG. 

图 1 给 出 了 系数 函数 C(s) 以 及 它 的 通 量 s. C(s) 的 图 形 . 当 C(s) 取 正 值 时 发 生前 向 扩散 ， 
用 来 实现 平滑 ， 当 C(s) 取 负 值 时 ， 发 生 后 向 扩散 ， 用 来 实现 边缘 的 锐 化 和 增强 . 


图 1: 关于 梯度 模 s 的 扩散 系数 函数 C(s) 及 通 量 s.C(s) 曲线 
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2.2 小 波 变换 
设 %(z, 切 是 二 维基 本 小 波 ， 则 任意 f(x,y) € L? (IR?) 的 连续 小 波 变 换 (CWT) 定义 为 加 


Wafle) = ff fehl dudo = (Fluo), bayt), (5) 
其 中 


s 表 示 尺 度 . 如果 池 (8, 切 是 可 分 离 的 ， 那 么 CWT 有 三 个 分 量 ， 即 水 平 细节 、 垂 直 细节 与 对 角 
细节 ， 本 文 主要 用 小 波 变换 提取 图 像 的 边缘 信息 ， 因 此 只 用 到 水 平 与 垂直 两 个 方向 的 分 量 ， 以 
下 讨论 中 ， 我 们 记 这 两 个 方向 的 小 波 变换 为 
WE f(z, y) = (f(u, v), YË „y(u, v)), k = 1, 2, 
EP yix, y) Sv? (cy) PH EAE SEAT LAR. 小 波 变换 的 模 定义 为 
M, f(x,y) = VIWI @ DP + WPF (ey). (6) 
车 取 尺度 s 为 {2};en， 则 称 相应 的 {WV 寺 f(z,9), W3 Fev) }jen X f(e, y) 的 二 进 小 波 变换 ， 二 
进 小 波 变换 的 优点 是 可 通过 卷 积 进行 快速 计算 . 
通常 用 CWT 检 测 图 像 边缘 ， 但 其 计算 比 二 进 小 波 变换 复杂 ， 因 此 本 文 在 以 下 扩散 模型 的 
构造 中 采用 二 进 小 波 变换 
图 2 给 出 了 应 用 小 波 变换 与 高 斯 梯度 的 含 噪 Lena 图 像 的 边缘 检测 效果 .图 2(a) 为 原始 图 
像 ， 图 2(b) 为 加 噪 图 像 ， 图 2(c) 和 图 2(d) 分 别 为 了 = 2,3 的 二 进 小 波 变换 模 图 像 ， 图 2(e) 和 
图 2(f) 分 别 为 经 过 方差 为 1 与 2 的 高 斯 函数 平滑 后 的 梯度 图 像 ， 由 图 可 见 ， 在 如 肩膀 与 镜面 的 
平滑 区 域内 ， 梯 度 图 像 比 小 波 模 图 像 含有 更 多 的 噪声 ， 理 论 上 ， 梯 度 图 像 的 噪声 随 高 斯 函数 的 
方差 增 大 逐渐 减 小 ， 但 同时 细小 边缘 也 逐渐 消失 . 


(a) (b) (c) 
(a) (e) (£) 


图 2: 小波 变 换 模 与 梯度 模 图 像 
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3 ”基于 小 波 变 换 模 的 FAB 扩 散 方 程 


3.1 基于 小 波 变换 模 的 FAB 扩 散 方 程 
根据 2.2 节 的 讨论 ， 利 用 小 波 变换 可 充分 提取 图 像 的 局 部 结构 信息 .因此 我 们 定义 如 下 的 基 
于 小 波 变换 的 FAB 扩散 模型 


| atent = div[C(Mz f(a,y))VF], te(0,+00), (2,y) ER2， 
f(x, y,0) = fo(x,y), 


其 中 C(s) 由 (2) 式 定义 , My f(x,y) 由 (6) 式 定义 , 该 模型 中 各 参数 的 作用 和 限定 条 件 同 FAB 模 
型 ， 与 Gilboa 的 FAB 扩散 模型 相 比 ， 扩 散 模 型 (7) 利用 小 波 变换 横 代 替 梯 度 模 确定 扩散 系数 ， 
因此 它 是 Gilboa 扩散 模型 的 推广 以 下 我 们 称 (7) 为 小 波 FAB 扩散 模型 ， 简 称 为 WFAB 扩散 
模型 . 

WFAB 扩散 模型 的 性 质 与 小 波 尺度 参数 7 有关、 当 尺 度 ; 较 小 时 ， 小 波 变 换 模 含 有 更 多 的 
小 边缘 信息 ， 同 时 也 含有 较 多 的 噪声 ， 扩 散 滤 波 对 图 像 边 缘 有 较 强 的 保持 和 增强 作用 ， 但 平滑 
能 力 降低 . 当 尺度 增 大 时 ， 小 波 变换 模 主要 包含 大 边缘 信息 ， 这 时 扩散 滤波 主要 增强 大 的 边 
缘 ， 对 噪声 和 细小 边缘 有 较 强 的 平滑 作用 . 因此 在 实际 应 用 时 应 适当 的 选择 尺度 参数 . 

3.2 基于 WEFAB 模型 的 扩散 算法 步骤 

步骤 1 选择 合适 的 小 波 和 尺度 参数 22， 对 f(z,y) = folz, y) 分别 做 沿 z 和 yy 方向 的 二 进 小 
波 变换 ， 得 WJ f(x,y),，W2 f(z,Yy); 

步骤 2 按 (6) 式 计算 小 波 变 换 模 Mas f(x,y): 

步骤 3 选择 满足 (3) HBR ky, ky, w, a, nn 和 ms; 

步骤 4 将 Mo; f(z,y) 和 参数 kp, ky, w, a, n, m, RA (2) 式 ， 计算 C(Mz; fle, y)): 

步骤 5 比较 前 -后 向 扩散 的 最 大 通 量 ， 如 果 


max{s-C(s)} > _ max {s-C(s)}, 











(7) 


ky -w<s<kytw 


DRO; 否则 转 步 又 3; 
BRE KE C(Mai fle, y) RA (7) A, 求解 该 定 解 问题 . 
在 用 和 迭代 格式 求解 (7) 时 ，C(s) 要 根据 前 一 次 滤波 结果 重新 计算 . 


4 数值 实验 


为 了 检验 WFAB 方法 的 平滑 增强 效果 ， 我 们 给 出 了 一 系列 数值 实验 结果 . 原始 图 像 是 大 小 
为 256 x 256 的 Boat 图 像 和 314 x 314 的 Tiffany 图像 ， 如 图 3(a) 和 图 4(a) 所 示 . 为 便于 比较 ， 
各 种 模型 的 求解 均 用 前 后 向 差分 混合 离散 格式 .实验 中 的 小 波 函 数 采用 二 次 B 样 条 小 波 . 

表 1 和 表 2 给 出 了 用 FAB 方 法 和 WFAB 方法 分 别 对 Boat Al Tiffany 图 像 进行 滤波 所 得 到 的 
滤波 图 像 的 峰值 信 噪 比 . 表 中 第 1 行为 输入 图 像 噪声 水 平 ， 其 中 ao 为 高 斯 白 噪声 的 标准 方 
#, input-PSNR 为 含 噪 图 像 的 峰值 信 噪 比 PSNR(dB)， 其 余 各 行为 输出 图 像 PSNR(dB). 在 
不 同 噪声 水 平 下 ，WEAB 方 法 的 去 噪 效果 均 优 于 传统 的 FAB 方法. 

图 3 和 图 4 分 别 给 出 了 用 WFAB 方法 和 FAB 方法 作 图 像 扩 散 滤波 的 滤波 图 像 ， 在 图 3 中 ， 
FAB 模型 的 参数 为 : ky = 30, ky = 180, w = 2, a = 0.08; 而 WFAB 模型 的 参数 为 : kj = 
30，k = 40, w = 5, a=0.3， 其 中 加 噪 图 像 的 噪声 水 平 为 c = 10. 图 4 中 ，FAB 模型 的 参 
BA: ky = 30, ky = 120, w = 2, a = 0.12，WFAB 模型 的 参数 为 : ky = 40, ky = 45, w= 
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2, a = 0.4， 其 中 加 噪 图 像 的 噪声 水 平 为 c = 20. 以 上 实验 的 时 间 步 长 均 取 为 0.1， 和 迭代 次 数 
为 5 次 . 





图 3: 前 向 -后 向 扩散 滤波 例子 ，(a) 原始 图 像 ; (b) SaNa RERE, REES Æ o = 10; 
(c) FAB 方 法 ; (d) WFAB 方 法， 其 中 j=3 





(a) (b) (c) (d) 


图 4: 前 向 -后 向 扩散 滤波 例子 ，(a) 原始 图 像 ; (b) 含 高 斯 白 噪声 图 像 ， 骂 声 标 准 方差 go = 20; 
(c) FAB 方法 ; (d) WFAB 方 法， 其 中 j=3 


表 1: Boat 图 像 的 FAB 方法 和 本 文 WFAB 方法 比较 








o/input-PSNR 

















| 10/32.04 20/29.33 30/28.48 40/28.09 50/27.88 
FAB 方 法 | 33.10 33.01 31.74 29.78 28.68 
WFAB 方 法 | 34.67 33.29 31.77 30.58 29.66 








R2: Tiffany 图像 的 FAB 方 法 和 本 文 WFAB 方法 比较 


o /input-PSNR 
10/32.09 20/29.39 30/28.53 40/28.19 50/27.95 


FAB 方法 | 33.99 33.89 32.33 29.94 28.78 
WFAB 方法 | 35.69 34.34 32.75 31.52 30.59 
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比较 而 言 ， 用 WFAB 方法 处 理 后 得 到 的 恢复 图 像 具 有 更 清楚 的 边缘 ， 如 Boat 图 像 中 的 船 
航 、 术 杆 以 及 天 空 的 白云 部 分 ， 边 缘 保 持 和 增强 效果 明显 . 这 说 明 ，WFAB 方法 在 处 理 边 缘 方 
面具 有 优势 ， 但 在 平滑 区 域 ， 如 Tiffany 手背 部 分 ， 不 如 FAB 方法 柔和 . 

为 了 进一步 说 明 WFAB 方法 的 增强 效果 ， 我 们 提取 了 Boat 和 Tiffany 恢复 图 像 上 的 一 段 局 
部 横 截 灰 度 值 曲线 ， 如 图 5 所 示 . 对 应 WFAB 方法 的 灰 度 曲线 在 边缘 处 明显 高 于 用 FAB 方法 
的 得 到 的 曲线 ， 这 说 明 WFAB 方法 对 边缘 的 增强 作用 强 于 FAB 方法 . 


























Ad: 恢复 图 像 的 横 截 灰 度 曲 线 ，(a) Boat 图 像 的 横 规 灰 度 值 曲 线 ; (b) Tiffany 图 像 的 
横 堆 灰 度 值 曲 线 ， 上 庶 线 对 应 于 FAB 方法 ， 实 线 对 应 WFAB 方法 


5 ”结论 


本 文通 过 将 小 波 变换 边缘 检测 的 思想 引入 到 FAB 扩 散 模 型 中 ， 分 析 研 究 了 基于 小 波 变换 
的 FAB 扩 散 模型 . 由 于 小 波 的 多 分 辨 性 ， 增 强 了 FAB 扩散 模型 的 抗 噪 能 力 和 边缘 增强 能 
同时 小 波 的 引入 也 为 扩散 系数 的 构造 提供 了 新 的 方法 .但 用 小 波 变换 检测 图 像 边 缘 时 易 产 生 边 
缘 偏 移 ， 从 而 引起 滤波 图 像 在 边缘 处 的 振动 ， 如 何 准确 定位 边缘 使 扩散 系数 更 准确 反映 图 像 的 
局 部 结构 信息 是 需要 进一步 探讨 的 问题 . 
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Forward-and-backward Diffusion Processes Based on Wavelet 
Transform Modulus 


ZHOU Feng, LIU Feng 


(Department of Information Sciences, School of Science, Xi’an Jiaotong University, Xi’an 710049) 


Abstract: This paper introduces the idea of the multiscale edge detection via wavelet modulus into 
the forward-and-backward (FAB) diffusion models. By utilizing the characteristics that the wavelet 
transform can well extract the local structure information of an image, we present a new method to 
construct an diffusivity using wavelet modulus at several scales, and establish the wavelet-based FAB 
diffusion model. Finally, we present several numerical experiments for image smoothing and enhancing 
to compare the proposed model with the previous FAB models. 
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